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regime based on a plurality of parameters, extracting pre- 
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and a fault level mapping for each of a plurality of fault 
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tion of the measurement vector onto the average direction of 
each of the plurality of fault types, determining a fault type 
based on which projection is maximum, and mapping the 
projection to a continuous-valued fault level using a lookup 
table. 
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METHOD AND APPARATUS FOR IN-SITU 
DETECTION AND ISOLATION OF 
AIRCRAFT ENGINE FAULTS 

STATEMENT REGARDING FEDERALLY 5 

SPONSORED RESEARCH & DEVELOPMENT 

This invention was made with Government support under 
contract number NAS 3-01135 awarded by the National 
Aeronautics and Space Administration (NASA). The Gov- to 
emment has certain rights in the invention. 

BACKGROUND OF THE INVENTION 

This invention relates to detecting and classifying faults in 15 
an operating machine, and more particularly to detecting and 
classifying faults in an operating aircraft engine using an 
Extended Kalman Filter architecture. 

Aircraft engines must maintain the highest achievable 
levels of reliability, because of their extreme safety-critical 20 
nature and because the vehicles powered by these engines 
represent enormous investments in resources. However, as 
with all machinery, small component failures and other 
operating faults may occur, owing to material failures, 
environmental disturbances, and normal deterioration dur- 25 
ing the operating life of an aircraft engine. 

Having faults go undetected and without compensating 
control actions can risk further damage and may accelerate 
deterioration, leading to higher safety risks. Similarly, when 
engine faults are detected by imprecise means and with high 30 
levels of uncertainty, operators are often obliged to take the 
most conservative measures, which typically involve abort- 
ing a takeoff or shutting down an engine during flight. Since 
these measures in themselves pose some risk to the aircraft 
and its occupants, it is important to be able to distinguish 35 
small faults for which more timely and less extreme mea- 
sures can safely be taken. 

Current engine health monitoring schemes detect only 
large faults and failures of the sensors, actuators, and control 
hardware. The architecture of the engine controls is based 40 
either in dual-redundant or tri-redundant hardware. Much of 
the diagnostic logic depends on comparing the redundant 
sensors to each other or to simple static models of the sensor. 
There is no systematic procedure for taking into account the 
behavior of the overall system by using a system model in 45 
concert with all of the available sensors. This causes current 
methods to be unable to detect faults until they reach a 
relatively large magnitude. Current monitoring also detects 
undesired and potentially damaging engine events like stalls 
and surges, but does not try to isolate the cause of the event. 50 

SUMMARY OF THE INVENTION 

An exemplary embodiment of the present invention 
includes a method for performing fault estimation based on 55 
residuals of detected signals, the method including compar- 
ing the detected signals with estimates of the detected 
signals, based on an extended Kalman filter, and outputting 
the residuals; and determining a fault type and a fault level 
by performing hypothesis testing on the residuals. 60 

In another exemplary embodiment of the present inven- 
tion, there is a method for performing a fault estimation 
based on residuals of detected signals, including: determin- 
ing an operating regime based on a plurality of parameters; 
extracting predetermined noise standard deviations of the 65 
residuals corresponding to the operating regime and scaling 
the residuals; calculating a magnitude of a measurement 


2 

vector of the scaled residuals and comparing the magnitude 
to a decision threshold value; extracting a mean direction 
and a fault level mapping for each of a plurality of fault 
types, based on the operating regime; calculating a projec- 
tion of the measurement vector onto the mean direction of 
each of the plurality of fault types; determining a fault type 
based on which projection is maximum; and mapping the 
projection to a continuous-valued fault level using a lookup 
table. 

In an additional exemplary embodiment of the present 
invention, there is a computer program product for enabling 
a computer to implement operations for performing a fault 
estimation based on residuals of detected signals, the com- 
puter program product comprising a computer readable 
medium and instructions on the computer readable medium, 
the operations including: determining an operating regime 
based on a plurality of parameters; extracting predetermined 
noise standard deviations of the residuals corresponding to 
the operating regime and scaling the residuals; calculating a 
magnitude of a measurement vector of the scaled residuals 
and comparing the magnitude to a decision threshold value; 
extracting a mean direction and a fault level mapping for 
each of a plurality of fault types, based on the operating 
regime; calculating a projection of the measurement vector 
onto the mean direction of each of the plurality of fault 
types; determining a fault type based on which projection is 
maximum; and mapping the projection to a continuous- 
valued fault level using a lookup table. 

In another exemplary embodiment of the present inven- 
tion, there is a method for detecting and isolating faults in a 
system, including: detecting a plurality of signals; determin- 
ing a residual of each of the plurality of signals; determining 
an operating regime based on a plurality of parameters; 
extracting predetermined noise standard deviations of the 
residuals corresponding to the operating regime and scaling 
the residuals; calculating a magnitude of a measurement 
vector of the scaled residuals and comparing the magnitude 
to a decision threshold value; extracting a mean direction 
and a fault level mapping for each of a plurality of fault 
types, based on the operating regime; calculating a projec- 
tion of the measurement vector onto the mean direction of 
each of the plurality of fault types; determining a fault type 
based on which projection is maximum; and mapping the 
projection to a continuous-valued fault level using a lookup 
table. In a further exemplary embodiment of the present 
invention, there is a computer program product for enabling 
a computer to implement operations for detecting and iso- 
lating faults in a system based on residuals of detected 
signals, the computer program product comprising a com- 
puter readable medium and instructions on the computer 
readable medium, the operations including: detecting a 
plurality of signals; determining a residual of each of the 
plurality of signals; determining an operating regime based 
on a plurality of parameters; extracting predetermined noise 
standard deviations of the residuals corresponding to the 
operating regime and scaling the residuals; calculating a 
magnitude of a measurement vector of the scaled residuals 
and comparing the magnitude to a decision threshold value; 
extracting a mean direction and a fault level mapping for 
each of a plurality of fault types, based on the operating 
regime; calculating a projection of the measurement vector 
onto the mean direction of each of the plurality of fault 
types; determining a fault type based on which projection is 
maximum; and mapping the projection to a continuous- 
valued fault level using a lookup table. 

In an additional exemplary embodiment of the present 
invention, there is an apparatus for detecting and isolating 
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faults in a system based on residuals of detected signals, the 
apparatus including: a processor configured to determine an 
operating regime based on a plurality of parameters; extract 
predetermined noise standard deviations of the residuals 
corresponding to the operating regime and scale the residu- 
als; calculate a magnitude of a measurement vector of the 
scaled residuals and compare the magnitude to a decision 
threshold: value; extract a mean direction and a fault level 
mapping for each of a plurality of fault types, based on the 
operating regime; calculate a projection of the measurement 
vector onto the mean direction of each of the plurality of 
fault types; determine a fault type based on which projection 
is maximum; and map the projection to a continuous- valued 
fault level using a lookup table. 

In another exemplary embodiment of the present inven- 
tion is a system for detecting and isolating faults based on 
residuals of detected signals, the system including: a detec- 
tor which detects the detected signals; an extended Kalman 
filter which compares the detected signals with estimates of 
the detected signals and outputs a plurality of residuals; and 
a processor which performs hypothesis testing on the residu- 
als to determine a fault type and a fault level. 

In another exemplary embodiment of the present inven- 
tion, there is a system for detecting and isolating faults based 
on residuals of detected signals, the system including: a 
detector which detects the detected signals; an extended 
Kalman filter which compares the detected signals with 
estimates of the detected signals and outputs a plurality of 
residuals; and a processor configured to determine an oper- 
ating regime based on a plurality of parameters; extract 
predetermined noise standard deviations of the residuals 
corresponding to the operating regime and scale the residu- 
als; calculate a magnitude of a measurement vector of the 
scaled residuals and compare the magnitude to a decision 
threshold value; extract a mean direction and a fault level 
mapping for each of a plurality of fault types, based on the 
operating regime; calculate a projection of the measurement 
vector onto the mean direction of each of the plurality of 
fault types; determine a fault type based on which projection 
is maximum; and map the projection to a continuous-valued 
fault level using a lookup table. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The advantages, nature and various additional features of 
the invention will appear more fully upon consideration of 
the illustrative embodiments of the invention which are 
schematically set forth in the figures, in which: 

FIG. 1 is a diagrammatical representation of a method for 
performing fault estimation based on residuals of detected 
signals according to an embodiment of the present inven- 
tion; 

FIG. 2 is a diagrammatical representation of a method for 
performing fault estimation based on residuals of detected 
signals according to another embodiment of the present 
invention; 

FIG. 3 is a diagrammatical representation of a method for 
detecting and isolating faults in a system according to an 
embodiment of the present invention; 

FIG. 4 is a diagrammatical representation of a system for 
detecting and isolating faults based on residuals of detected 
signals, according to an embodiment of the present inven- 
tion; and 

FIG. 5 is a diagrammatical representation of data pro- 
vided by an operating regime lookup table. 
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DETAILED DESCRIPTION OF THE 
INVENTION 

The present invention will be explained in further detail 
5 by making reference to the accompanying drawings, which 
do not limit the scope of the invention in any way. 

Modem aircraft engines employ full-authority digital con- 
trols, which make use of sensors deployed throughout the 
engine. This invention describes how these same sensor 
to measurements can be used to monitor the health of the 
engine, including its actuators and the sensors themselves. 
By using a system model with available sensors, this inven- 
tion is able to isolate the cause of the engine event when it 
occurs. Moreover, this invention is able to distinguish small 
15 faults from large faults. Along with the safety advantages of 
being able to distinguish small faults from larger faults, the 
ability to detect small faults early enables more timely 
engine maintenance, which reduces costs and extends the 
operating life of the engine. 

20 The invention will now be taught using various exemplary 
embodiments. Although the embodiments are described in 
detail, it will be appreciated that the invention is not limited 
to just these embodiments, but has a scope that is signifi- 
cantly broader. The appended claims should be consulted to 
25 determine the true scope of the invention. Prior to describing 
the embodiments in detail, however, the meaning of certain 
terms will be explained. 

One embodiment of this invention resides in a computer 
system. Here, the term “computer system” is to be under- 
30 stood to include at least a memory and a processor. In 
general, the memory will store, at one time or another, at 
least portions of an executable program code, and the 
processor will execute one or more of the instructions 
included in that executable program code. It will be appre- 
35 ciated that the terms “executable program code,” “software,” 
and “instructions” mean substantially the same tiling for the 
purposes of this description. It is not necessary to the 
practice of this invention that the memory and the processor 
be physically located in the same place. That is to say, it is 
40 foreseen that the processor and the memory might be in 
different physical pieces of equipment or even in geographi- 
cally distinct locations. 

The above-identified invention may be embodied in a 
computer program product, as will now be explained. On a 
45 practical level, the software that enables the computer sys- 
tem to perform the operations described in detail further 
below may be supplied on any of a variety of media. 
Furthermore, the actual implementation of the approach and 
operations of the invention may actually be statements in a 
50 computer language. Such computer language statements, 
when executed by a computer, cause the computer to act in 
accordance with the particular content of the statements. 
Furthermore, the software that enables a computer system to 
act in accordance with the invention may be provided in any 
55 number of forms including, but not limited to, original 
source code, assembly code, object code, machine language, 
compressed or encrypted versions of the foregoing, and any 
and all equivalents now known or hereafter developed. 

One familiar with this field will appreciate that “media”, 
60 or “computer-readable media”, as used here, may include a 
diskette, a tape, a compact disc, an integrated circuit, a 
ROM, a CD/DVD, a cartridge, a memory stick or card, a 
remote transmission via a communications circuit, or any 
other medium useable by computers, including those now 
65 known or hereafter developed. For example, to supply 
software for enabling a computer system to operate in 
accordance with the invention, the supplier might provide a 
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disc or might transmit the software in some form via satellite 
transmission, via a direct wired or a wireless link, or via the 
Internet. Thus, the term, “computer readable medium” is 
intended to include all of the foregoing and any other 
medium by which software may be provided to a processor. 5 

Although the enabling software/code/instructions might 
be “written on” a disc, “embodied in” an integrated circuit, 
or “carried over” a communications circuit, it will be appre- 
ciated that, for the purposes of this discussion, the software 
will be referred to simply as being “on” the computer to 
readable medium. Thus, the tenn “on” is intended to encom- 
pass the above mentioned and all equivalent and possible 
ways in which software can be associated with a computer 
readable medium. 

For the sake of simplicity, therefore, the term “program 15 
product” is thus used to refer to a computer readable 
medium, as defined above, which has on it any form of 
software to enable a computer system to operate according 
to any embodiment of the invention. 

Having explained the meaning of various terms, the 20 
invention will now be described in detail, in the context of 
a method. 

In an exemplary embodiment of the invention, there is a 
method for performing fault estimation based on residuals of 
detected signals. As illustrated in FIG. 1, the method 25 
includes: comparing the detected signals with estimates of 
the detected signals, based on an extended Kalman filter 
(EKF), and outputting the residuals (step 100 ); and deter- 
mining a fault type and a fault level by performing hypoth- 
esis testing on the residuals (step 101 ). Examples of the 30 
detected signals correspond to actual sensor measurements. 
The EKF compares actual sensor measurements to estimates 
provided by an internal model and outputs error signals, i.e., 
residuals. The EKF is described in the following references: 
Athans, M. (1996), The Control Handbook , pp. 589-594, 35 
CRC Press, United States and Anderson, B. D. O., Moore, 

J. B., Optimal Filtering , Prentice-Hall, Englewood Cliffs 
N.J., 1979. 

Prior to operation of the exemplary embodiments of the 
method for performing fault estimation, a training process is 40 
implemented in which an estimator is trained offline. An 
engine model is used in the training process to determine the 
noise variances and the final values of the residuals for each 
of the representative fault types and levels, including the 
no-fault case. To account for variation in engine behavior 45 
over the flight envelope, which may be defined by ambient 
temperature, altitude, mach number, and thrust level, the 
domain of possible variation in the flight envelope param- 
eters is divided into representative regimes. The final values 
of the sensor residuals are logged for each of these regimes 50 
by averaging over all test cases pertaining to a particular 
regime. The sensor noise variance is similarly segregated by 
regime. However, the noise is assumed to be independent of 
the faults. Thus, the computation of the standard deviations 
is performed only on the no-fault training data, after segre- 55 
gating the data by regime. 

The training process enables computing the noise stan- 
dard deviations, the mean directions for faults, the mappings 
between fault levels and vector magnitudes, and the assign- 
ment of a decision threshold value. In an exemplary embodi- 60 
ment, the fault level mappings are stored in a lookup table, 
based upon the operating regimes. The training process is 
repeated for each sensor, for each fault type, and for each 
regime comprising the flight envelope. 

In an exemplary embodiment of the invention, the fore- 65 
going method is implemented in an aircraft to detect and 
isolate errors in the aircraft engine and its associated actua- 
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tors and sensors. For an engine operating normally, the 
residuals are small and contain only measurement noise. 
When a fault occurs, the residuals respond in a maimer 
systematic to the type and severity of the fault. The type of 
fault and the level of the fault are determined by performing 
hypothesis testing on the residuals. For example, a Bayesian 
Hypothesis Test may be performed, in which the likelihoods 
of various predefined fault types are assessed given the 
current residuals. A decision on the fault type is made, which 
is followed by a correlation computation to determine the 
fault level or severity. 

In the Bayesian Hypothesis Testing, the various fault 
types and levels, as well as the no-fault condition, are 
expressed in terms of their signatures in a measurement 
space defined by the scaled EKF residuals, where the scaling 
is performed by dividing each residual by the standard 
deviation of that signal’s noise level. This space is of 
dimension p, where p is the number of sensors. The values 
attained by the residuals when the various fault conditions 
are imposed are thus represented as positions in the space, 
and these positions are compared to real-time measurements 
to asses the health of the engine. Bayesian Hypothesis 
Testing is described in the following reference: Van Trees, 
H. L., Detection, Estimation, and Modulation Theory, John 
Wiley & Sons, N.Y., 1968 (Sections 2.1 thru 2.4). 

Because the sensors are subjected to noise and unknown 
biases, the fault conditions give rise not to discrete positions 
in the measurement space, but rather to probabilistic distri- 
butions in the space. More precisely, these are conditional 
probability functions defined on the p-dimensional space, 
given the various fault hypotheses. If the sensor noise is 
considered to be Gaussian and white, the fault hypotheses 
may be represented as ellipsoidal functions centered on 
various mean positions in the space, with the axes of the 
ellipsoid sized according to the noise variance in each 
dimension. When the noise is independent across sensors, it 
is customary to normalize the dimensions of the measure- 
ment space by the standard deviation of the corresponding 
sensor noise, so that the ellipsoidal probability density 
functions degenerate to spherical functions of uniform 
radius. Further, in the absence of a priori knowledge of 
engine faults, the various fault hypotheses are assumed to be 
equally likely over a given time interval. Under these 
assumptions, the implementation of an optimal Bayesian 
Hypothesis test that minimizes the probability of error (false 
positives, false negatives, and misclassifications) is accom- 
plished by means of a simple distance computation between 
the real-time residuals and the p-dimensional reference 
hypotheses, after dividing each residual by the standard 
deviation of the sensor noise. This is because the conditional 
probabilities are represented directly by distance, in the 
standard Euclidian sense, within the normalized measure- 
ment space. 

In another exemplary embodiment of the present inven- 
tion, there is a computer program product for enabling a 
computer to implement the operations for performing fault 
estimation based on residuals of detected signals, the com- 
puter program product comprising a computer readable 
medium and instructions on the computer readable medium, 
the operations including: comparing the detected signals 
with estimates of the detected signals, based on an extended 
Kalman filter (EKF), and outputting the residuals (step 100); 
and determining a fault type and a fault level by performing 
hypothesis testing on the residuals (step 101), as described 
above in relation to FIG. 1. 

FIG. 2 illustrates an exemplary embodiment of the present 
invention, in which there is a method for performing a fault 
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estimation based on residuals of detected signals. This 
method includes: determining an operating regime from a 
plurality of parameters (step 201); extracting predetermined 
noise standard deviations of the residuals corresponding to 
the operating regime and scaling the residuals (step 202); 5 
calculating a magnitude of a measurement vector of the 
scaled residuals (step 203); determining if the magnitude is 
at or above the decision threshold value (step 204); extract- 
ing a mean direction and a fault level mapping for each of 
a plurality of fault types, based on the operating regime (step to 
205); calculating a projection of the measurement vector 
onto the mean direction of each of the plurality of fault types 
(step 206); determining a fault type based on which projec- 
tion is maximum (step 207); and mapping the projection to 
a continuous-valued fault level using a lookup table (step 15 
208). An example of data that can be provided in an 
operating regime lookup table is illustrated in FIG. 5, which 
depicts the magnitude of a projection versus severity of 
fault. However, the operating regime lookup table and the 
data therein are not limited to the illustration in FIG. 5. 20 

The detected signals may include any or all of actuator 
signals, sensor signals and engine signals. 

In step 201, an operating regime is determined from a 
plurality of parameters. 

In step 202, predetermined noise standard deviations of 25 
the residuals are extracted from a source of data, i.e., a 
lookup table, for example, and the residuals are scaled. In an 
exemplary embodiment, the scaling includes normalizing 
the residuals. 

In step 203, the magnitude of the measurement vector of 
the residuals is calculated. The measurement vector may be 
determined by dividing each residual by a noise standard 
deviation. This magnitude is compared to the decision 
threshold value, which is a predetermined value. The selec- 
tion of the decision threshold value is based upon a balance 
between the false alarm rate and the fault detection rate. Too 
low a decision threshold value increases the false alarm rate, 
while too high a decision threshold value increases the 
likelihood that actual faults will go undetected. 

40 

In step 204, it is determined whether the magnitude is at 
or above the decision threshold value. If it is determined that 
the magnitude is below the decision threshold value, it is 
determined that there is no fault. The operating regime 
represents one of a plurality of possible portions of the 45 
region in which a signal may be present. The parameters 
used to determine the operating regime may include flight 
envelope parameters. For example, the flight envelope 
parameters may include, but are not limited to ambient 
temperature, altitude, mach number, and thrust level. 50 

After it is determined that there is a fault of some type, a 
mean direction and a fault level mapping are extracted for 
each of the plurality of fault types based on the operating 
regime, in step 205. The fault types may include, but are not 
limited to, one or more of a sensor fault, an actuator fault, 55 
a first machine fault, and a second machine fault. The mean 
directions are unit vectors approximating the contours 
defined by the end values of the residuals, which can be used 
in a correlation computation to determine the fault type. In 
an exemplary embodiment, the mean direction for each fault go 
type includes a set of p-dimensional vectors, where p 
represents a number of sensors. The fault level mapping may 
be determined via a lookup table that associates fault levels 
with length along the appropriate fault contour. 

In step 206, a calculation of a projection of the measure- 65 
ment vector onto the mean direction of each of the plurality 
of fault types is performed. Based on a maximum projection 
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of the measurement vector onto the mean direction of each 
of the plurality of fault types (step 207), the fault type is 
determined. 

For the determined fault type, the projection is mapped to 
a continuous-valued fault level, using a lookup table (step 
208). The continuous values may be obtained by interpo- 
lating or extrapolating from predetermined fault levels. 

FIG. 3 illustrates a method for detecting and isolating 
faults in a system. The method illustrated in FIG. 3 corre- 
sponds to the method illustrated in FIG. 2, but further 
includes the steps of detecting signals (step 301) and deter- 
mining the residuals of the detected signals (step 302). Once 
steps 301 and 302 are performed, the method of FIG. 3 
follows the method illustrated in FIG. 2. Since the method 
of FIG. 2 is described above, the description of these steps 
is not repeated here. 

In another exemplary embodiment of the present inven- 
tion, there is a computer program product for enabling a 
computer to implement operations for detecting and isolat- 
ing faults in a system based on residuals of detected signals, 
the computer program product including a computer read- 
able medium and instructions on the computer readable 
medium, the operations including: detecting signals (step 

301) ; determining residuals of the detected signals (step 

302) ; determining an operating regime from a plurality of 
parameters (step 201); extracting predetermined noise stan- 
dard deviations of the residuals corresponding to the oper- 
ating regime and scaling the residuals (step 202); calculating 
a magnitude of a measurement vector of the scaled residuals 
(step 203); determining if the magnitude is at or above the 
threshold value (step 204); extracting a mean direction and 
a fault level mapping for each of a plurality of fault types, 
based on the operating regime (step 205); calculating a 
projection of the measurement vector onto the mean direc- 
tion of each of the plurality of fault types (step 206); 
determining a fault type based on which projection is 
maximum (step 207); and mapping the projection to a 
continuous-valued fault level using a lookup table (step 
208). Since these steps are described above, the description 
is not repeated here. 

FIG. 4 illustrates a block diagram of a system for detect- 
ing and isolating faults based on residuals of detected 
signals, according to an embodiment of the present inven- 
tion. As shown in FIG. 4, the system 400 includes a detector 
401 which detects signals; an extended Kalman filter 402 
which compares the detected signals with estimates of the 
detected signals and outputs a plurality of residuals; and a 
processor 403 which performs hypothesis testing on the 
residuals to determine a fault type and a fault level. The 
processor 403 may be configured to operate in real time, i.e., 
during the operation of the system, without introducing a 
delay into the operation of the system. Also, the hypothesis 
testing may be Bayesian Hypothesis Testing. 

The detector 401 may include a plurality of sensors 404, 
which are disposed in predetermined locations throughout 
the system 400. The sensors 404 are configured to monitor 
a machine, which, in the present embodiment, is an aircraft 
engine 405. In a particular embodiment, sensors 404 are 
disposed on the machine. Sensors 404 may also be config- 
ured to monitor the actuators 406 or other subsystems within 
the system 400. In particular embodiments, sensors 404 are 
disposed on the actuators or other subsystems. The engine 
controller 407 of FIG. 4 provides control signals to actuators 
406, which control operations of the engine 405. 

In an exemplary embodiment of the present invention, the 
processor 403 is configured to calculate a magnitude of a 
measurement vector of the residuals and compare the mag- 
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nitude to a decision threshold value; determine an operating 
regime based on a plurality of parameters; extract a mean 
direction and a fault level mapping for each of a plurality of 
fault types, based on the operating regime; calculate a 
projection of the measurement vector onto the mean direc- 
tion of each of the plurality of fault types; determine a fault 
type based on which projection is maximum; and map the 
projection to a continuous-valued fault level using a lookup 
table. These operations of the processor are described above 
in relation to FIG. 2. 

While the invention has been described in terms of 
specific embodiments, those skilled in the art will recognize 
that the invention can be practiced with modification within 
the spirit and scope of the claims. Namely, although the 
present invention has been discussed in the context of 
aircraft engine applications, it is contemplated that the 
present invention can be employed in all applications in 
which faults of a machine are detected and classified. 

What is claimed is: 

1. A method for performing a fault estimation based on a 
plurality of residuals of a plurality of detected signals from 
a machine, comprising: 

determining an operating regime based on a plurality of 
parameters; 

extracting predetermined noise standard deviations of 
said residuals corresponding to said operating regime 
and scaling said residuals; 

calculating a magnitude of a measurement vector of said 
scaled residuals and comparing said magnitude to a 
decision threshold value; 

extracting a mean direction and a fault level mapping for 
each of a plurality of fault types, based on said oper- 
ating regime; 

calculating a protection of said measurement vector onto 
said mean direction of each of said plurality of fault 
types; 

determining a fault type based on a maximum projection 
of said measurement vector onto said mean direction of 
each of said plurality of fault types; and 

mapping said projection to a continuous- valued fault level 
using a lookup table, wherein said calculating the 
magnitude of the measurement vector and comparing 
the same to the decision threshold value, determining 
the fault type based on the maximum projection of the 
measurement vector, and mapping said projection to 
the continuous-valued fault level are provided for sub- 
sequent use in perfonning said fault estimation of the 
machine. 

2. The method of claim 1, wherein, if said magnitude is 
below said decision threshold value, it is determined that 
there is no fault. 

3. The method of claim 1, wherein said parameters 
comprise flight envelope parameters. 

4. The method of claim 1, wherein said scaling comprises 
normalizing said residuals by dividing said residuals by said 
noise standard deviations. 

5. The method of claim 1, further comprising correlating 
said measurement vector with said mean direction of each of 
said plurality of fault types. 

6. The method of claim 1, wherein said residuals comprise 
extended Kalman filter residuals. 

7. The method of claim 1, wherein said detected signals 
comprise at least one of actuator signals, sensor signals, and 
engine signals. 

8. The method of claim 1, wherein said detected signals 
comprise actuator signals, sensor signals, and engine sig- 
nals. 
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9. The method of claim 1, wherein said fault types 
comprise at least one of a sensor fault, an actuator fault, a 
first machine fault and a second machine fault. 

10. The method of claim 1, wherein said fault types 
5 comprise a sensor fault, an actuator fault, a first machine 

fault and a second machine fault. 

11. The method of claim 1, wherein said mean direction 
for each fault type comprises a set of p-dimensional vectors, 
where p represents a number of sensors. 

to 12. A computer program product for enabling a computer 
to implement operations for perfonning a fault estimation 
based on a plurality of residuals of a plurality of detected 
signals from a machine, the computer program product 
comprising a computer readable medium and instructions on 
15 the computer readable medium, the operations comprising: 

detennining an operating regime based on a plurality of 
parameters; 

extracting predetennined noise standard deviations of 
said residuals corresponding to said operating regime 
20 and scaling said residuals; 

calculating a magnitude of a measurement vector of said 
scaled residuals and comparing said magnitude to a 
decision threshold value; 

extracting a mean direction and a fault level mapping for 
25 each of a plurality of fault types, based on said oper- 
ating regime; 

calculating a projection of said measurement vector onto 
said mean direction of each of said plurality of fault 
types; 

30 determining a fault type based on a maximum projection 
of said measurement vector onto said mean direction of 
each of said plurality of fault types; and 

mapping said projection to a continuous-valued fault level 
using a lookup table, wherein said calculating the 
35 magnitude of the measurement vector and comparing 
the same to the decision threshold value, determining 
the fault type based on the maximum projection of the 
measurement vector, and mapping said projection to 
the continuous-valued fault level are provided for sub- 
40 sequent use in perfonning said fault estimation of the 
machine. 

13. The computer program product of claim 12, wherein, 
if said magnitude is below said decision threshold value, it 
is detemiined that there is no fault. 

45 14. The computer program product of claim 12, wherein 

said parameters comprise flight envelope parameters, 
wherein said detected signals comprise at least one of 
actuator signals, sensor signals, and engine signals, and said 
system comprises an engine, wherein said fault types com- 
50 prise at least one of a sensor fault, an actuator fault, a first 
machine fault and a second machine fault. 

15. The computer program product of claim 12, wherein 
said scaling comprises normalizing said residuals by divid- 
ing said residuals by said noise standard deviations, and 

55 wherein said residuals comprise extended Kalman filter 
residuals. 

16. The computer program product of claim 12, wherein 
the operations further comprise: 

detecting a plurality of signals; and 
60 detennining a residual of each of said plurality of signals, 
and 

wherein the fault estimation is performed for a system that 
includes an aircraft engine. 

17. A method for detecting and isolating faults in a 
65 system, comprising: 

detecting a plurality of signals; 

detennining a residual of each of said plurality of signals; 
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determining an operating regime based on a plurality of 
parameters; 

extracting predetermined noise standard deviations of 
said residuals corresponding to said operating regime 
and scaling said residuals; 5 

calculating a magnitude of a measurement vector of said 
scaled residuals and comparing said magnitude to a 
decision threshold value; 

extracting a mean direction and a fault level mapping for 
each of a plurality of fault types, based on said oper- 1° 
ating regime; 

calculating a projection of said measurement vector onto 
said mean direction of each of said plurality of fault 
types, 

determining a fault type based on a maximum projection 15 
of said measurement vector onto said mean direction of 
each of said plurality of fault types; and 

mapping said projection to a continuous-valued fault level 
using a lookup table wherein said calculating the mag- 
nitude of the measurement vector and comparing the 20 
same to the decision threshold value, determining the 
fault type based on the maximum projection of the 
measurement vector, and mapping said projection to 
the continuous-valued fault level are provided for sub- 
sequent use in perfonning said fault estimation of the 25 
machine. 

18. The method of claim 17, further comprising correlat- 
ing said measurement vector with said mean direction of 
each of said plurality of fault types; 

wherein said measurement vector is determined by divid- J ' 
ing each residual by a noise standard deviation; 

wherein said residuals comprise extended Kalman filter 
residuals; and 

wherein, if said magnitude is below said decision thresh- 35 
old value, it is determined that there is no fault. 

19. The method of claim 17, wherein said parameters 
comprise flight envelope parameters, wherein said detected 
signals comprise at least one of actuator signals, sensor 
signals, and engine signals, and said system comprises an 40 
engine, wherein said fault types comprise at least one of a 
sensor fault, an actuator fault, a first machine fault and a 
second machine fault, and wherein said system comprises an 
aircraft engine. 

20. An apparatus for detecting and isolating faults in a 45 
system based on a plurality of residuals of a plurality of 
detected signals from the system, said apparatus comprising: 

a processor configured to detennine an operating regime 
based on a plurality of parameters, extract predeter- 
mined noise standard deviations of said residuals cor- 50 
responding to said operating regime and scale said 
residuals, calculate a magnitude of a measurement 
vector of said scaled residuals and compare said mag- 
nitude to a decision threshold value, extract a mean 
direction and a fault level mapping for each of a 55 
plurality of fault types, based on said operating regime, 
calculate a projector of said measurement vector onto 
said mean direction of each of said plurality of fault 
types, detennine a fault type based on a maximum 
projection of said measurement vector onto said mean 60 
direction of each of said plurality of fault types and map 
said projection to a continuous-valued fault level using 
a lookup table, wherein the processor is configured to 
calculate the magnitude of the measurement vector and 
compare the same to the decision threshold value, 65 
determine the fault type based on the maximum pro- 
jection of the measurement vector, and map said pro- 
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jection to the continuous-valued fault level for subse- 
quent use in detecting and isolating the faults in the 
system. 

21. The apparatus of claim 20, wherein said processor is 
configured to operate in real time. 

22. The apparatus of claim 20, wherein said apparatus is 
disposed in an aircraft and said system comprises an aircraft 
engine. 

23. The apparatus of claim 20, wherein said processor 
comprises an aircraft engine controller. 

24. The apparatus of claim 20, wherein said apparatus is 
configured to detect and isolate faults in a plurality of 
actuators, a plurality of sensors, and an engine. 

25. A system for detecting and isolating faults based on a 
plurality of residuals of a plurality of detected signals, said 
system comprising: 

a detector which detects said detected signals; 

an extended Kalman filter which compares said detected 
signals with estimates of said detected signals and 
outputs a plurality of residuals; and 

a processor which performs hypothesis testing on said 
residuals to determine a fault type and a fault level, 
wherein said processor is configured to determine an 
operating regime based on a plurality of parameters, 
extract predetermined noise standard deviations of said 
residuals corresponding to said operating regime and 
scale said residuals, calculate a magnitude of a mea- 
surement vector of said scaled residuals and compare 
said magnitude to a decision threshold value, extract a 
mean direction and a fault level mapping for each of a 
plurality of fault types, based on said operating regime, 
calculate a projection of said measurement vetcor onto 
said mean direction of each of said plurality of fault 
types, detennine a fault type based on a maximum 
projection of said measurement vector onto said mean 
direction of each of said plurality of fault types, and 
map said projection to a continuous-valued fault level 
using a lookup table. 

26. The system of claim 25, wherein said detector com- 
prises a plurality of sensors. 

27. The system of claim 25, wherein said hypothesis 
testing comprises Bayesian hypothesis testing. 

28. The system of claim 25, wherein said processor is 
configured to operate in real time. 

29. The system of claim 25, wherein said processor 
comprises an aircraft engine controller and is disposed in an 
aircraft and said system comprises an aircraft engine, and 
wherein said system is configured to detect and isolate faults 
in a plurality of actuators, a plurality of sensors, and an 
engine. 

30. A method for perfonning fault estimation based on a 
plurality of residuals of a plurality of detected signals from 
a machine, said method comprising: 

comparing said detected signals with estimates of said 
detected signals, based on an extended Kalman filter, 
and outputting said residuals; and 

detennining a fault type and a fault level by perfonning 
hypothesis testing on said residuals, said detennining 
further comprises determining an operating regime 
based on a plurality of parameters, extracting predeter- 
mined noise standard deviations of said residuals cor- 
responding to said operating regime and scale said 
residuals, calculating a magnitude of a measurement 
vector of said scaled residuals and compare said mag- 
nitude to a decision threshold value, extracting a mean 
direction and a fault level mapping for each of a 
plurality of fault types, based on said operating regime, 
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calculating a projection of said measurement vetcor 
onto said mean direction of each of said plurality of 
fault types, determining a fault type based on a maxi- 
mum projection of said measurement vector onto said 
mean direction of each of said plurality of fault types, 5 
and mapping said projection to a continuous -valued 
fault level using a lookup table, wherein said calculat- 
ing the magnitude of the measurement vector and 
comparing the same to the decision threshold value, 
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determining the fault type based on the maximum 
projection of the measurement vector, and mapping 
said projection to the continuous-valued fault level are 
provided for subsequent use in performing said fault 
estimation of the machine. 

31. The method of claim 30, wherein said hypothesis 
testing comprises Bayesian hypothesis testing. 



